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Abstract 

Background: Canine distemper virus (CDV) infects a variety of carnivores, including wild and domestic Canidae. In 
this study, we sequenced and phylogenetic analyses of the hemagglutinin (H) genes from eight canine distemper 
virus (CDV) isolates obtained from seven raccoon dogs {Nyctereutes procyonoides) and a giant panda {Ailuropodo 
melanoleuca) in China. 

Results: Phylogenetic analysis of the partial hemagglutinin gene sequences showed close clustering for geographic 
lineages, clearly distinct from vaccine strains and other wild-type foreign CDV strains, all the CDV strains were 
characterized as Asia-1 genotype and were highly similar to each other (91.5-99.8% nt and 94.4-99.8% aa). The giant 
panda and raccoon dogs all were 549Y on the HA protein in this study, irrespective of the host species. 

Conclusions: These findings enhance our knowledge of the genetic characteristics of Chinese CDV isolates, and 
may facilitate the development of effective strategies for monitoring and controlling CDV for wild canids and 
non-cainds in China. 
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Background 

CD is caused by the canine distemper virus (CDV) which 
belongs to the Morbillivirus genus of the Paramyxoviridae 
virus family. This virus infects a broad range of animals, 
such as Mustelidae (ferrets, minks, skunks, weasels, and 
badgers), Procyonidae (raccoons), Ursidae (bears and 
pandas), Viverridae (civets, genets, and linsangs), hyaenidae 
(hyenas), and Felidae (lions and tigers) [1-7]. The gen- 
ome of CDV encodes the following genes: matrix (M), 
fusion (F), hemagglutinin (H), nucleocapsid (N), poly- 
merase (L), and phosphoprotein (P). The H gene protein 
is responsible for viral attachment to the cell host [8] 
and is the most variable protein described for all members 
of the genus Morbillivirus [9]. Experimental investigation 
showed that changing one amino acid (Y549H) in the 
CDV-H protein of one dog strain can decreased expression 
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of specialist traits and increased expression of generalist 
traits, thereby confirming its functional importance [10]. 

The H gene has been used to investigate the genetic 
relationships among the various strains [11-13]. On the 
basis of the full-length sequence of the H gene of CDV 
strains identified globally, at least seven main geographic 
groups (genotypes) appear to circulate among the vari- 
ous susceptible hosts, namely America- 1, America-2, 
Asia-1 and Asia-2, Europe, Europe wildlife and Arctic 
[14-17]. In addition, some CDV strains identified in Africa, 
Asia and Argentina [18-20] appear to diverge substantially 
and might represent separate geographic groups. 

N glycosylation has been shown to be important for 
the correct folding, transport, and function of other 
paramyxovirus fusion and attachment glycoproteins [21]. 
There are significant molecular differences between the 
glycoproteins of wild-type and vaccine CDV strains. It has 
been hypothesized that reduced N glycosylation is an im- 
portant attenuating factor and that an increase in N glyco- 
sylation may eventually result in vaccine failure [15,22], 
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Although the use of live attenuated vaccines has 
greatly helped to prevent CDV infections in susceptible 
animals, the dental abnormalities and enamel erosion 
that are found in many young and adult giant pandas 
[23,24] may be related to early infection with CDV 
[25,26]. Several episodes also have been reported in rac- 
coon dogs [20]. So in this work we analyzed a total of 7 
samples from raccoon dogs with clinical signs of canine 
distemper and 1 sample from giant panda which were 
submitted to our laboratories and the viruses were ana- 
lyzed using a fragment of the H gene. 

Results 

CDV detection by NP-based RT-PCR 

We analyzed a total number of 108 samples taken 
from raccoon dogs with suspect of CDV infection. 
CDV RNA was detected in 56 out of 108 (51.8%) 
blood samples, most of the raccoon dogs, between the 
age of 15-30 weeks. Negative samples were assayed 
twice to confirm the results. Forty-one of the 56 posi- 
tive clinical specimens (73.2%) were obtained from 
raccoon dogs which had been vaccinated against CDV 
at least once. Twelve samples (21.4%) were obtained 
from unvaccinated raccoon dogs, and the vaccination 
status of the three remaining raccoon dogs was un- 
known. Forty-nine of the 56 positive animals (87.5%), 
showed symptoms typical of CDV, such as respiratory 
disease and neurological dysfunction, two raccoon 
dogs were asymptomatic, while no records were avail- 
able for the other five raccoon dogs. The giant panda 
showed a mild fever (anal temperature 39.0°C). 

Sequence and phylogenetic analysis of the H gene 

Amplification of fragments of the expected size was 
obtained with H gene and the sequences determined. 
After removing primer-derived sequences and some se- 
quences with inconsistencies at the 5' and 3' ends, a 
594 bp-long fragment of the H gene were obtained. 

The neighbor joining tree (Figure 1) generated using 
a 594 base pair (bp) long partial sequences of the H 
gene (amplified with primers "B")> evidenced geographic 
patterns (genotypes) resembling the patterns described 
previously using the full-length H gene (Asia 1, Asia 2, 
Europe, Arctic, Americal and America 2). In the phylo- 
genetic tree (Figure 1) all the Chinese wild-type viruses 
(8/8) grouped together in one branch in the Asia-1 
genotype. In this genotype, other Chinese CDV strains 
identified from domestic dogs or other wildlife species 
after the 1990s were also included, along with some 
Japanese and Korean strains detected before 1998. The 8 
Asia-1 strains displayed identity to each other in the H 
gene (91.5-99.8% nt and 94.4-99.8% aa), while identity to 
the vaccine strains CDV3, Onderstepoort and Lederle 
was lower (85.4-89.2% nt and 84.4-88.3% at aa level) 
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Figure 1 Phylogenetic relationships among 73 CDV isolates 
based on the amino acid sequences of the H gene. Distance 
values were calculated by the ClustalW program within the MEGA 
5.0 software package. The Neighbor-Joining algorithm was used to 
generate the tree. Bootstrap values were calculated on 1000 
replicates. Only bootstrap values >75% were shown. Filled 
quadrate (■) indicates the 8 Chinese wild-type CDV strains 
analyzed in this study. 



(Table 1). The 4 strains (Sichuan 02-05) displayed high 
identity to each other (97.4-99.8% nt and 98.4-99.8% aa). 
On the other hand, although strain Sichuan06 appeared 
to belong to the main Asia 1 clade, it is more distant 
from the group (Figure 1). 
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Table 1 Homology of nucleotide and amino acid sequences of CDV H gene 

SichuanOI Sichuan02 Sichuan03 Sichuan04 Sichuan05 Sichuan06 Sichuan07 Sichuan08 CDV3 Onderstepoot Lederle 

Nucleotides(%) 



SichuanOI 




93.1 


94.1 


94.9 


94.9 


94.6 


98.3 


98.0 


86.5 


85.4 


86.7 


OlLI lUdl lUZ 


QA A 




Q9. ^ 


Q9. 3 


Q9. 3 

70.J 


yQ.Q 


QA Q 
yH.y 


QA Q 


RQ 1 
oy. I 


9.7 1 


9.Q 1 
oy.Z 


Sichuan03 


95.0 


98.4 




97.4 


97.4 


94.6 


95.2 


92.1 


87.7 


86.5 


88.0 


Sichuan04 


95.0 


99.4 


99.7 




99.8 


97.3 


97.5 


91.7 


88.0 


86.5 


88.2 


Sichuan05 


95.0 


99.5 


99.6 


99.8 




97.5 


97.7 


91.9 


88.2 


86.6 


88.5 


Sichuan06 


95.5 


98.9 


99.4 


99.5 


99.4 




95.9 


91.5 


88.0 


86.5 


88.2 


Sichuan07 


99.8 


94.4 


95.0 


95.0 


95.0 


95.5 




93.9 


88.0 


86.4 


88.3 


Sichuan08 


99.7 


94.5 


95.3 


95.0 


95.0 


95.5 


99.5 




88.4 


87.2 


88.6 


CDV3 


86.0 


86.6 


87.2 


87.2 


87.5 


87.7 


86.0 


86.4 




96.1 


99.0 


Onderstepoot 


84.4 


84.9 


85.5 


85.5 


85.7 


86.0 


84.4 


84.7 


93.9 




96.0 


Lederle 


86.6 


87.2 


87.7 


87.7 


87.9 


88.3 


86.6 


86.8 


97.2 


93.9 





Amino acids(%) 



Analysis of the amino acid sequences of the H gene of 
wild-type CDV strains 

The inferred aa sequences of the partial H protein of 
the 8 wild-type Chinese CDV strains and of vaccine 
(America- 1 genotype) strains were aligned. The giant 
panda and raccoon dogs all were 549Y on the HA 
protein from the current study (Table 2). Inspection of 
the H protein alignment (Table 3) revealed that all the 
7 Asia-1 strains possessed 9 potential N-linked glycosyla- 
te sites (N-X-S/T) at aa positions 19-21, 149-151, 
309-311, 391-393,422-424, 456-458, 584-586, 587-589 
and 603-605. In the H protein of strain SichuanOI, the 
sites at positions 584-586 and 603-605 were disrupted. 
The strain Onderstepoort had only four glycosylation 
sites, while strains CDV3, Lederle and Convac displayed 
seven glycosylation sites. 

Discussion 

Despite the vaccination procedures adopted in China in 
the last decades, CDV is still a serious threat to breed- 
ing raccoon dogs as well as to domestic dogs and it is 



responsible for relevant economic losses in Chinese fur in- 
dustry [20] . In order to determine the frequency of this in- 
fection, we analyzed the samples taken from raccoon dogs 
with suspect of CDV infection. Our result revealed a rela- 
tively high (51.8%) prevalence of CDV infection animals 
sampled with clinical signs consistent with CDV infection. 
This indicated that CDV still represents a high risk to the 
canine population in China. Furthermore, all samples were 
tested twice and the results were reproducible. One pos- 
sible explanation for the 52 animals thought to be infected 
with CDV were found negative by PCR screening is the 
ability of detecting low copy numbers of RNA. Another 
possible explanation for the negative is the presence of Taq 
polymerase inhibitors on blood samples, since it is known 
that blood material contains nucleases that degrades RNA, 
and substances, such as heparin and hemoglobin were 
found to inhibit the activity of Taq DNA polymerase [27]. 

The molecular analysis of the amplifications obtained 
by the use of primer pair "B" helped to clarify the rela- 
tionship of the China strains with other CDV strains 
reported in other parts of the world. The analysis has 



Table 2 CDV isolates from China, organs from which they were isolated and their GenBank accession numbers 



Serial no 


Isolate 


Host 


Foster mode 


Vaccinated 


Organ source 


GenBank accession no. 


549 


1 


SichuanOI 


Raccoon dog 


Captivity 


Yes 


Blood 


JX406733 


Y 


2 


Sichuan02 


Raccoon dog 


Captivity 


Yes 


Blood 


JX406734 


Y 


3 


Sichuan03 


Raccoon dog 


Captivity 


Yes 


Blood 


JX406736 


Y 


4 


SichuanOI 


Raccoon dog 


Captivity 


Yes 


Blood 


JX406737 


Y 


5 


Sichuan05 


Raccoon dog 


Captivity 


Yes 


Blood 


JX406738 


Y 


6 


Sichuan06 


Raccoon dog 


Captivity 


Yes 


Blood 


JX406739 


Y 


7 


Sichuan07 


Raccoon dog 


Captivity 


Yes 


Blood 


JX406740 


Y 


8 


Sichuan08 


Giant Panda 


Captivity 


Yes 


Blood 


JX406735 


Y 
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Table 3 Distribution status of N-glycosylation sites of H gene in different genetypes of CDV 

Genetype CDV isolate N-glycosylation sites in each sequence 



19 149 309 391 421 456 584 587 603 



Asia-1 


SichuanOI 


N 


N 


N 


N 


N 


N 


N 


N 


N 




Sichuan02 


N 


N 


N 


N 


N 


N 


N 


N 


N 




Sichuan03 


N 


N 


N 


N 


N 


N 


N 


N 


N 




Sichuan04 


N 


N 


N 


N 


N 


N 


N 


N 


N 




Sichuan05 


N 


N 


N 


N 


N 


N 


N 


N 


N 




Sichuan06 


N 


N 


N 


N 


N 


N 


N 


N 


N 




Sichuan07 


N 


N 


N 


N 


N 


N 


N 


N 


N 




Sichuan08 


N 


N 


N 


N 


N 


N 


N 


N 


N 


Vaccines 


Onderstepoort 


N 


N 






N 






N 






Convac 


N 


N 




N 


N 


N 




N 


N 




CDV3 


N 


N 




N 


N 


N 




N 


N 




Lederle 


N 


N 




N 


N 


N 




N 


N 



Note:"N": N-glycosylation Sites;"-": negative. 



shown that the strains of giant panda and raccoon dogs in 
this study are placed on one branches of the phylogenetic 
tree, were characterized as Asia-1 genotype. The H gene- 
based neighbor- joining phylogenetic tree analysis using se- 
lected sequences retrieved from GenBank revealed 
geographic-related patterns of segregation. Our findings 
are consistent with previous studies demonstrating the oc- 
currence of strains belonging to genetically distinct CDV 
lineages in Asia [20,28,29]. The Chinese Asia-1 strains 
displayed a high genetic conservation with other Asia-1 
CDVs, which detected from raccoon dogs or other carni- 
vores in China over nearly two decades (Figure 1). The 
strain Sichuan08 was obtained from a mild fever giant 
panda in the Mar of 2012, and it showed the identity be- 
tween Sichuan02-06 were 94.5-95.5%. Because there were 
few reports of giant panda clinical onset of CDV infection 
up to now, so we couldn't determine whether or not the 
fever is the clinical manifestations of the giant panda. 

The observed diversity between the vaccine strains 
and recent wild-type CDVs may be accounted for by sev- 
eral mechanisms, such as adaptation to new host species 
[30,31], antigenic escape [32-35] and/or genetic recom- 
bination between wild- type strains [36], variously driving 
the evolution of the virus. It has been described that 
changes in N-glycosylation of the H protein may affect 
neutralization by antibodies and replication in vitro 
[15,37,38]. So, the connected aspartic amide N glycosyla- 
tion site potentially is a spotlight in H proteins between 
vaccine and wild strains of CDV. Usually, there are four 
(Onderstepoort strain) or seven (CDV3, Lederle and 
Convac strain) potential sites in the vaccine strains. 
However, seven or nine sites have been detected in all 
wild CDV strains, of which 309-311 N-connected amide 
asparagine glycosylation sites are specific to CDV field 



strains. Some studies believed that the variants from H 
protein glycosylation played a crucial role in the anti- 
genic differences and increase in N glycosylation may 
eventually result in vaccine failure [15,39]. So, we 
conjectured that the giant panda and raccoon dogs 
infected CDV maybe due to this reason in this time. 

Based on phylogenetic analyses of complete sequences 
data from CDV strains retrieved from domestic dog and 
non-dog species it was hypothesized that positive selec- 
tion drives residues 530 and 549 of the HA gene, and 
that the presence of specific residues at these positions 
resulted in the emergence of CDV as a disease of novel 
host species [31]. Recently, it has been suggested that 
there no evidence that amino acid 530 was strongly af- 
fected by host species, while wild canids there was a 
nearly significant trend towards a 549Y bias, non-canid 
strains showed no significant bias towards either H or Y 
at site 549 [9]. In this study, We found that amino acid 
549 Y was conserved within CDV lineages, regardless of 
host species. CDV transmission in wild carnivore and 
non-canids species may also most often occur between 
individuals within a species, and will be influenced by a 
range of factors such as population size and ranging pat- 
terns [40,41]. 

Conclusions 

In conclusion, the results presented by this study from a 
relatively small number of samples obtained from two host 
species conform to information obtained by other studies 
in China, that information currently available on CDV 
strains in China, including the results from the current 
study, is insufficient to determine whether or not the Asia 
1 genotype is the predominant genotype in China. Analysis 
of CDV strains from a variety of host species will provide a 
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more in-depth understanding of the global ecology of 
CDV, So, detect these more research on a broader range of 
potential carnivore hosts is required. 

Materials and methods 

Viruses and clinical specimens 
Ethics statement 

The animal from which specimens were collected, was 
handled in accordance with animal protection law of the 
People s Republic of China (a draft of an animal protec- 
tion law in China released on September 18, 2009). This 
study was approved by the National Institute of Animal 
Health Animal Care and Use Committee at Sichuan 
Agricultural University (approval number 2010-020). 

From 2011 to 2012, 108 samples (whole blood) of rac- 
coon dogs with suspect of CDV infection were sent to 
our laboratory for diagnostic confirmation and were 
identified only in 56 samples (51.8%) by using a RT- 
PCR, which amplifies a 242 bp-long fragment of the nu- 
cleoprotein (NP) gene with primers "A" ( Table 4). Out 
of 56 samples, 7 could be included in this study based 
on the amount of the sample and the amount of RNA 
extracts available for additional analysis. The samples 
were submitted to our laboratory from 3 different pri- 
vate fur farms within yaan, Sichuan province. The or- 
gans of the giant panda (whole blood) were submitted to 
our laboratory by China Conservation and Research 
Center for the Giant Panda, Yaan, China. Detailed infor- 
mation on the origin and the accession numbers of the 
CDV-positive samples is shown in Table 2. 

Extraction of nucleic acids 

Viral RNA was extracted from leucocytes obtained by 
washing with nuclease-free water. Briefly, 500 \A of 
nuclease-free water were added to 500 \A of whole blood 
in a 1.5 microcentrifuge tube, and centrifuged 7 min at 
5,000 g. A total of 500 [A of the supernatant was 
discarded, then 500 \A of nuclease-free water were added, 
repeating this step 3 times. Total RNA was extracted 
with TRI Reagent (Sigma®, USA), following the manufac- 
turer s instructions. 

RT-PCR 

Reverse transcription was carried out by PrimeScript™ 
RT reagent Kit (TaKaRa, Dalian, China) according to the 



manufacturers instruction: 37°C for 15 min, 85°C for 5 sec, 
using 7 \A of RNA, 2 (il 5 x PrimeScript buffer, 0.5 (il ran- 
dom primer and 0.5 \A PrimeScript RT Enzyme. The CDV- 
positive samples were tested with primers "B" (Table 4), 
which amplifying a 613 bp fragment (corresponding to po- 
sitions 8006-8619 of the Onderstepoort strain genome). 
The thermo-cycling condition were as follows: 5 min at 
94°C (initial denaturation), 30 cycles of 30 sec at 94°C, 
30 sec at 51°C, 30 sec at 72°C, and a final extension of 
10 min at 72°C. 

Cloning and sequencing 

PCR products of the correct size (613 bp in length) were 
amplified and cloned into the pMD 19-T vector (TaKaRa, 
Dalian, China). For each CDV strain, 3-5 positive recom- 
binant plasmids were sequenced in both directions using 
primer M13 (Shanghai Invitrogen Biotechnology Com- 
pany, Shanghai, China) and additional primers selected on 
the basis of the obtained sequences. The accession num- 
bers are shown in Table 2. 

Amino acid analysis 

To investigate the differences between the Chinese wild- 
type CDV strains and the vaccine strains available in the 
market in China, the deduced amino acid (aa) sequences 
of the H protein were aligned. Multiple sequence align- 
ment was carried out using the software package 
MegAlign (Lasergene7.0, DNAStar Inc., Madison, WI). 
The aa present at sites 549 on the HA protein was deter- 
mined for the eight strains from the current study. The 
potential N-linked glycosylation sites of the H protein 
were determined using the software NetNGlycl.O. 

Phylogenetic analysis 

The eight China CDV sequences were compared to those 
of 64 CDV strains that were obtained from various part of 
the world up until November 2010. The Accession num- 
bers for these H gene nucleotide sequences were obtained 
from the GenBank nucleotide database at the National 
Center for Biotechnology Information (NCBI). The nucleo- 
tide sequences were analyzed preliminary with cognate se- 
quences available in GenBank using BLAST software 
(http://www.ncbi.nlm.nih.gov/). Multiple sequence align- 
ment was made with ClustalW software [42] included in 
MEGA 5.0 software [43]. 



Table 4 Description of the primer pairs 



Primer pair 


Sequence 5'-3' 


Target 


Sense 


Purpose 


Amplicon size (bp) 


A 


CGGAAATCAACGGACCTAAAT 


N 


+ 


Diagnostic 


242 




TCOTGAGCmCGACCOT 










B 


TGGTOACAAGATGGTATO 


H 


+ 


Phylogenetic 


613 




CAACACCAC TAAATOGACT 











Guo et al. Virology Journal 2013, 10:109 
http://www.virologyj.eom/content/10/1/109 



Page 6 of 7 



Sequences for comparison were obtained corresponded 
to nucleotide sequences of the H gene. Phylogenetic ana- 
lysis of nucleotide were computed using the Kimura two 
parameter algorithm and the tree was constructed with the 
Neighbor Joining method with a 1000 bootstrap repeats 
using the MEGA 5.0 Software [43]. The measles virus 
(MeV) (NC001498), phocine (seal) distemper virus (PDV) 
(AF479274), peste-des- petits -ruminants virus (PPRV) 
(NC006383), hendra virus (HeV) (NC001906) and dolphin 
morbillivirus (DMV) (NC005283) were used as outgroup 
in the neighbor joining trees. 
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